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Abstract 

Coulomb energy differences between mirror nuclei with A < 7Q are cal¬ 
culated within the framework of the nuclear shell model using an effective 
Coulomb plus isovector and isotensor interaction. Absolute binding energies 
for proton-rich nuclei are predicted by adding the calculated Coulomb shifts 
to experimentally measured binding energies for the neutron-rich mirror. The 
location of the proton drip-line is investigated, as well as candidates for the 
exotic decay mode known as di-proton emission. Taking into account the 
lifetimes of competing decay modes and limits imposed by experimental se¬ 
tups, it is concluded that the best candidates for the observation of correlated 

di-proton emission are ^^Fe, '^^Ni, and ®^Se. 
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I. INTRODUCTION 


The structure of exotic nuclei, i.e., nuclei with extreme isospin values, is one of the 
most exciting challenges in low-energy nuclear physics today. Detailed theoretical studies 
of exotic nuclei, when confronted with experiments, will yield important information about 
the interaction between nucleons in the nucleus and the validity of our models for the 
structure of nuclei. In addition, the study of exotic nuclei is essential to many fundamental 
issues in physics today; in particular, the weak interaction and nuclear astrophysics. For 
example, it is believed that many of the heavy elements in the universe are produced by the 
radiative capture of neutrons (r-process) [Q] or protons (rp-process) p| on unstable nuclei. 
The competition between beta decay and particle capture traces out a path that synthesizes 
the known elements. The details of this path, and, hence, the abundance of the elements 
produced, depends on the temperature of the site, as well as explicit nuclear properties such 
as binding energies, level densities, spectroscopic factors, and beta-decay lifetimes. 

An additional new feature of proton-rich nuclei that will be explored in the next few years 
is the possibility of a new decay mode known as diproton emission. Because of the pairing 
interaction, a nucleus with an even number of protons {Z, N) is generally more tightly bound 
than a. {Z — 1,N) nucleus, but, because of the symmetry energy and Coulomb repulsion, 
it may be unbound relative to the {Z — 2, N) system. The number of candidates for the 
observation of this decay mode, however, is sharply limited by the two-proton separation 
energy. This is in part due to the fact that emission is a competing decay mechanism, 
and, because of the large Coulomb energy difference between the parent and daughter nuclei, 
the beta-decay lifetimes are of the order 1-100 ms. In addition, a further constraint on the 
observation of diproton emission can be imposed by the experimental apparatus, since, in 
many experiments, the parent nucleus must live long enough to be identified. Generally 
speaking, these two practical constraints limit the observable lifetime for diproton decay to 
10“® — 10“^ s. On the other hand, the decay rate for diproton emission is determined by 
the probability to penetrate through the Coulomb barrier, which, in turn, is exponentially 
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dependent on the two-proton separation energy. As will be shown here, the number of 
candidates for which the observation of diproton decay is practical, is limited to nuclei with 
two-proton separation energies between 0.9 and 1.4 MeV. 

One of the principal motivations for the construction of radioactive beam facilities is to 
study the properties of nuclei near the limits of stability. Very few nuclei near the proton 
drip-line have been identihed, and the heaviest and most proton-rich nucleus observed to 
date is ^®Ni Q. Even more difficult than the identihcation of an exotic nucleus is the 
measurement of its mass, and, at present, predictions regarding the r- and rp-process, must 
rely on theoretical estimates for nuclear binding energies. 

Several methods have been used to obtain theoretical estimates for absolute binding 
energies. One is the liquid-drop formula and associated variants, such as the microscopic- 
macroscopic approach . In general, these models are determined by htting a set of liquid- 
drop parameters while including effects due to pairing and shell corrections to experimental 
data over a wide range of nuclei, and have been found to reproduce known nuclear masses 
at the level of approximately 800 keV (see for example Ref. [^). Although the microscopic- 
macroscopic approach gives a good global description of nuclear binding energies and is the 
method of choice for heavy nuclei where detailed microscopic calculations are not feasible, 
there are notable discrepancies between experimental and calculated binding energies with 
neutron number between 20 and 40, as is illustrated in Fig. 1 of Ref. |^. 

For lighter nuclei, however, more accurate binding energies can be achieved using the 
nuclear shell model, since, in many cases, it is necessary only to compute the Coulomb 
energy difference between mirror nuclei p-|^. In this paper. Coulomb energy differences 
are computed for mirror nuclei in the fp shell for 46 < A < 70. By then making use 
of experimental data for the neutron-rich members tabulated in Ref. |p, absolute binding 
energies are predicted with an estimated accuracy at the level of 50 — 200 keV. With these 
binding energies, two-proton separation energies are computed and rough estimates for the 
lifetimes for diproton decay are made. Given the practical constraints on the decay halflives 
mentioned above, it is found that the best candidates for the experimental observation of 
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correlated two-proton emission are '^^Fe, ^®Ni, and ®^Se. 

This paper is organized into five sections. In Section y, the systematics of Coulomb 
energy differences between analog nuclei are discussed, while in Section H a shell-model 
description of these energy shifts is presented. Candidates for the exotic decay mode known 
as di-proton emission are presented and analyzed in Section and concluding remarks are 
collected in Section 0. 


II. SYSTEMATICS OF COULOMB ENERGY DIFFERENCES 


If the nuclear Hamiltonian is composed of only one- and two-body parts, quite generally, 
it may be separated into three components. The dominant part, which is also responsible 
for most of the nuclear binding energy, is due to the strong interaction and is isoscalar in 
nature. The other two components are due to both the Coulomb interaction and charge 
non-symmetric parts of the nucleon-nucleon interaction, and are isovector and isotensor in 
character. If the isovector and isotensor components are weak relative to the isoscalar com¬ 
ponent, then the binding energies for the members of an isospin multiplet may be obtained 
within the context of the isobaric mass multiplet equation (IMME) 

BE{A, T, T„ t) = a{A, T, t) + b{A, T, z)T, + c(A, T, z)T|, (1) 


where T and = {Z — N)/2 denote the isospin and its third component for the members 
of the isospin multiplet, Z, N, and A = Z + N are the number of protons, neutrons, and 
nucleons, respectively. The label i in Eq. (|l|) represents all other quantum numbers needed 
to denote the state, such as angular momentum, state number, etc. The coefficients a, b, 


and c separately depend on the isoscalar isovector, and isotensor components of the 
nuclear Hamiltonian, respectively. 

In shell-model calculations, the isoscalar part of the nuclear Hamiltonian is usually de¬ 
termined empirically by fitting to experimental binding energies and levels that have had 
the Coulomb energy subtracted off in an average way (c.f. |]I^-[1T[]). The predictive power 
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of these effective Hamiltonians is indicated by the rms deviation between experimental data 
and calculated binding energies, and, to date, the best empirically determined Hamiltonian 
is that due to Wildenthal for use in the 0 ^ 5 / 2 , 0 ^ 3 / 2 , and lsi /2 orbitals, where the rms 
deviation between theory and experiment is of the order 200 keV. For the most part, this 
interaction may be thought of as indicative of the best accuracy that may be achieved within 
the framework of the nuclear shell model. For other model spaces, such as the /p-shell (de- 
hned by the O/ 7 / 2 , O/ 5 / 2 , IP 3 / 2 , and lpi /2 orbitals), the effective shell-model Hamiltonian is 
less well determined and the deviation between theory and experiment is somewhat larger 
and is of the order 300 keV [0 . With this in mind, we must conclude that any attempt to 
compute absolute binding energies from hrst principles in the shell model would include an 
uncertainty of at least 200-300 keV in the a coefficient of Eq. ©• 

From Eq. (0), the binding energy difference between iosbaric analogs with = ±T is 
given by 


BE{A,T,T, = T,t)-BE{A,T,T, 


-T,i) = 2b{A,T,i)T. 


( 2 ) 


Therefore, the most accurate way to predict absolute binding energies for proton-rich nuclei 
whose analog has an experimentally measured mass is to compute the 6 -coefficient for the 
multiplet (or the Coulomb energy difference) and add 2bT to the experimental binding 
energy, BEf.xp{A,T,Tz = —T,i), of the neutron-rich analog. The overall uncertainty in the 
predicted binding energy is then of the order 

6BE{A, T, T, = T,t) = ^{26bTY + 6{BE,^p{A,T,T, = -T,i))\ (3) 


where 5b is the uncertainty in the 6 -coefhcient and 5BEexp{A,T,Tz = —T,i) is the uncer¬ 
tainty in the experimental binding energy. In many cases, it is possible to estimate the 
6 -coefhcient with an uncertainty of the order 30-40 keV |^, and, therefore, it may be possi¬ 
ble to predict the binding energies of extreme proton-rich nuclei at the level of 100-200 keV. 

In Refs. |§-|§], the procedure outlined above was used to predict the absolute binding 
energies of nuclei with 36 < H < 55. Those three works approached Eq. (0) using slightly 
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different methods, but with about the same level of accuracy, as is indicated by the overall 
agreement between them. In Ref. P, Eq. (|^) was evaluated using shell-model calculations 
for the pure O/ 7 / 2 -shell nuclei and a weak coupling approximation for those nuclei that 
spanned both the 0 ^ 3/2 and O/ 7/2 orbits. In Ref. 0, all the Coulomb energy differences 
were evaluated within the framework of the shell model using the 0 ^ 3/2 and O/ 7/2 orbits 
and an empirical isospin-nonconserving (INC) interaction [^. Finally, in Ref. [§], Eq. 
was evaluated using a method based on a parameterization of the Coulomb displacement 
energies [|I^ . The overall success of these works, and agreement between them, is essentially 
due to their empirical foundations. In each, a set of parameters was £t to experimental 
data, and the models were then extrapolated to predict the masses of unknown nuclei. This 
work is an extension of Ref. [0 in which absolute binding energies of proton-rich nuclei with 
46 < A < 70 are predicted by computing the Coulomb displacements within the framework 
of the nuclear shell model. 

Before continuing with the details and the results of the shell-model calculations, it is in¬ 
structive to examine the systematic behavior of the Coulomb displacement energies. Indeed, 
one of the reasons for the success of the three different methods is the smooth behavior as 
a function of nucleon number A exhibited by experimental ^-coefficients. In addition, for a 
given mass number, the fe-coefficients are essentially constant to within 100 keV or so, as can 


be seen from Tables 3-7 in Ref. |]^. This behavior is easily understood from the liquid-drop 
model, where the Coulomb energy of a sphere of radius R = with charge Ze is given 

by 

3 {Zef 


Ec = 


5 R 


The Coulomb energy difference between analog nuclei is then 

3 3e^ , 9 / 0 , 

=- 2AT = -A^/^(2T), 

5 R 5 ro 


( 4 ) 


( 5 ) 

( 6 ) 


where A = Z + N and the isospin is defined by T = \Z — N\/2. Hence, by comparing 
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Eqs. (H) and (||), it is seen that the fe-coefficient is expected to increase as . 

For a comparison with experimental ^-coefficients, we tnrn to the more sophisticated 
liqnid-drop parameterization of Ref. P|. Here the form of the “Conlomb” energy will be 
ontlined, and the Conlomb energy difference between analog nnclei will be evalnated nsing 
the parameters dehned in Ref. 0. In macroscopic models, the “Coulomb” contribution to 
the binding energy is 

Ecoul = - C 4 ^ + fikfTp)— - Ca{N - Z), (7) 


where Ci = 3e^/5ro and C 4 = 5/4(3/27r)^/^ci. The hrst two terms in Eq. are the direct 
and exchange Coulomb energies, the third is the proton-form factor correction, and the last 
is the charge-asymmetry energy. The factor B 3 is the shape-dependent relative Coulomb 
energy, which, to leading order for a spherical shape, is given by 


5 75 105 


( 8 ) 


with yo = {ro/aden)A^^^ rs-/ 1.657y4^/^. The proton form factor / is dependent on the Fermi 
wave number kf = (97rZ/4H)^/^(l/ro) and the proton rms radius = 0.8 fm (see Eq. (8) 
of Ref. [Q), and for nuclei with H ~ 50 and Z ~ A/2 may be accurately approximated by 
/ = —0.214 MeV. Using tq = 1.16 fm and Ca = 0.145 MeV from Ref. |Q, the 6 -coefficient 
derived from the Coulomb-energy difference between analog nuclei is 


6 ^^ = [0.7448H2/3 _ 1.882 + 1.535H-^/^ - 0.782871"^] MeV. (9) 


Shown in Fig. |l| is a comparison between Eq. (^) (solid line) and experimental 6 - 
coefficients (solid squares). The experimental data comprise 116 6 -coefficients, and were 
taken from Tables 3-7 in Ref. and the known ground-state analog mass differences tabu¬ 
lated in Ref. ^ for A < 59. In the hgure, all the 6 -coefficients for a given mass number were 
averaged together, and error bar reflects both the standard deviation and the experimental 
uncertainties. Generally speaking, for a given A, the 6 -coefficients are roughly constant, with 
the mean standard deviation being 61 keV. From the microscopic point of view, deviations 
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from constancy can be expected for two reasons. First, in some cases single-particle orbits 
from different major oscillator shells come into play, as in 15 < A < 17, and second, near the 
limits of stability, the single-particle orbits are nearly nnbonnd, and the Thomas-Erhman 
shift [^,^] needs to be acconnted for (see Section [TTTI) . 

From Fig. |^, it is evident that the experimental ^-coefficients exhibit a global be¬ 
havior. On the other hand, Eq. (P) tends to nnderestimate the fe-coefficients for A < 40, and 
the rms deviation with the data is 138 keV. Within the context of a global parameterization, 
a slight improvement on Eq. (P) can be obtained by fitting to the experimental data, and 
an rms deviation of 102 keV is achieved with 

b = [0.710712/3 - 0.946] MeV, (10) 

which is also represented in Fig. |I| by the dashed line. For the most part, Eq. o leads to 
global description of the Conlomb energy differences between analog nnclei with an accnracy 
of the order 100|Z —A^| keV. To improve npon this, it is necessary to acconnt for local nnclear 
strnctnre via a microscopic model, which is the topic of the next section. 


III. SHELL-MODEL CALCULATIONS OF COULOMB ENERGY DIFFERENCES 


In this section, the procednre for compnting the Conlomb energy difference between ana¬ 
log nnclei within the framework of the shell model is ontlined. In Ref. ||^ , empirical isovector 
and isotensor, or isospin-nonconserving (INC), Hamiltonians were determined for several 
shell-model spaces by constraining them to reprodnce experimental b- and c-coefficients. 
The primary components of the empirical interactions were the Coulomb interaction and 
two-body isovector and isotensor interactions. In general, the empirical two-body isovector 
interaction was rather weak, while the isotensor interaction was found to to be consistent 
with the differences observed in the proton-proton and proton-neutron scattering lengths. 
The deviations between theoretical and experimental b and c coefficients were of the order 
of 30 keV and 15 keV, respectively. 
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In this work, proton-rich nuclei in the mass range 46 < A < 70 are investigated. For 
all but two cases (the T = 1/2, A = 69 isodoublet and the T = 1, A = 70 isotriplet), the 
binding energy of the neutron-rich analog has been measured and is tabulated in Ref. [Q. 
The shell-model calculations were performed using the shell-model code OXBASH in 
proton-neutron formalism using the conhguration space dehned by the O/ 7 / 2 , O/ 5 / 2 , IP 3 / 2 , 
and lpi /2 orbitals (the fp shell) and the FPD 6 Hamiltonian given in Ref. [|T^. Here, instead 
of computing the 6 -coefhcients, the “Coulomb” energy differences between analogs were 
computed directly by adding the INC interaction to the FPD 6 Hamiltonian. Due to the 
large dimensions present, some truncations on the model space were found to be necessary. 
For A < 59, all conhgurations contained within the O/ 7/2 and IP 3/2 orbits were included, 
and the truncations were based on the number of particles permitted to be excited out of 
the 0 / 7/2 and IP 3/2 orbits into the O/ 5/2 and lpi /2 orbits. Generally, this ranged from two to 
four particles so that the total dimensions (with good angular momentum) were less than 
14,000. For A > 60, the O/ 7/2 orbit was taken to be a closed core, and the O/ 5 / 2 , IP 3 / 2 , and 
lpi /2 single-particle energies of the FPD 6 interaction were modihed so as to reproduce the 
levels of ®^Ni under this assumption. For the most part, it was found that the Coulomb 
energy shifts were not particularly sensitive to the applied truncations, and in a few cases 
where the effects of the truncations were tested, differences of only a few keV were found. 
Therefore, the applied model-space truncations are not expected to provide a signihcant 
contribution to the uncertainty in the computed binding energies. 

The INC interaction used here consists of a Coulomb plus charge-dependent (isotensor) 


interaction, and is described in detail in Ref. |^. An important parameter for this interac¬ 
tion is the oscillator frequency, hu, since the Coulomb components are scaled as a function 
of A by the factor 0 


^(A) = 


huj{A) 

11.096 


1/2 


( 11 ) 


Generally, hui is chosen to reproduce experimental rms charge radii, and for many nuclei it 
can be accurately parameterized by 
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huj{A) = 45A-^/3 _ 25A-2/3 MeV. 


(12) 


It is important to note, however, that for A > 45, Eq. (p!2D underestimates huj as compared to 
values derived from experimental charge radii. Indeed, in Ref. the value of 10.222 MeV 
was used for A = 53 as opposed to the value of 10.208 MeV implied by Eq. (^). In 
addition, the /p-shell INC interaction was reht in Ref. 0 . where it was found that better 
overall agreement between theoretical and experimental Coulomb energy shifts was obtained 
using oscillator frequencies derived from the rms charge radii of Hartree-Fock calculations 
using the Skyrme M* interaction. These values of huo are tabulated in Table |, and are used 
in the present work. 

In Refs. [^,^, the htted INC interactions were able to reproduce the experimental 
6 -coefficients for /p-shell nuclei with an rms deviation of approximately 33 keV. However, 
the most difficult parameters to determine for the INC interaction are the Coulomb single¬ 
particle energies for the O/ 5/2 and lpi /2 orbits, as there is very little experimental data 
available that is sensitive to these quantities. In Ref. 0 . these single-particle energies were 
£t upon by making assumptions regarding spin assignments for excited levels in ®^Cu and 
®®Zn. In retrospect, these levels are probably not appropriate for determining parameters 
for heavier nuclei because of uncertainties in spin assignments and the fact that the lev¬ 
els comprising the assumed doublet at 1.040 MeV in ®^Cu are unbound, and strong 
Thomas-Ehrman shifts may apply (see below). Also, shell-model calculations for 

the J'^ = 1/2“ and 5/2“ states in ®®Zn indicate that these levels are predominantly lp |/2 
conhgurations. On the other hand, the the beta-endpoint energies for both ®^Ga and 
®®As are sensitive to the O/5/2 Coulomb single-particle energy and were used to help £x 
this parameter. In regards to the lpi /2 single-particle energy, however, no data exists that 


will definitively set this parameter. For this reason, the value obtained in Ref. which 
also happens to reproduce the 6 -coefficients for the assumed 1/2“ states in A = 57 and 59, 
is used. 

Two additional concerns that affect this work are: (1) whether the /p-shell alone is 
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sufficient to describe the nuclei in question, and (2) the effect of the Thomas-Ehrman shift on 
the Coulomb displacement energies near the drip line. As a measure of the appropriateness 
of just the /p-shell for the calculations, we examine the the excitation energies of the hrst 
J = 2+ states m. N = Z, even-even nuclei in the region 60 < A < 80. Shown in Table 
are the experimental PB| , PB| excitation energies of the these states in comparison with the 
values obtained with the FPD6* interaction (FPD6 modihed as indicated above after closing 
the 0 / 7/2 orbit). Overall, there is good agreement between the calculated and experimental 
values until A = 76, where there is a sudden drop in the excitation energy, which is an 
indication of the onset of collective behavior that would necessitate the inclusion of orbits 
from the next major shell, such as the Qgg /2 orbit. Given the results in Table |I|, the /p-shell 
is sufficient to describe the nuclei studied in this work. 

In general. Coulomb energies are computed using harmonic oscillator, or sometimes 
bound Woods-Saxon, single-particle wave functions for the protons, with the length scale 
chosen to reproduce experimental rms charge radii. Near the drip line, however, this ap¬ 
proximation can be inadequate. Because they are loosely bound, the proton single-particle 
wavefunctions are pushed out of the nuclear interior, and, as a consequence, the Coulomb 
energy is reduced. This shift in the Coulomb energy was hrst noted by Thomas [M] and 


Ehrman in the A = 13 system, and is most important for light nuclei where the Coulomb 
barrier, which acts to conhne the wave function in the nuclear interior, is smaller, and for 
orbits with little or no centrifugal barrier, eg., the Si /2 orbitals. This effect is well illustrated 
by the single-particle states in A = 17, where the Coulomb displacement energy of the 
= 5/2+ (the Qd ^/2 orbit) ground state is 3.543 MeV, while the shift for the = 1/2+ 
state, which is a lsi /2 single-particle state that is bound by only 107 keV, is 3.168 MeV. 
The inhuence of the centrifugal barrier is also apparent in these nuclei, as the Coulomb shift 
for the J"" = 3/2 state (the 0 ^ 3/2 spin-orbit partner of the ground state), which is unbound 
by 4.5 MeV, is 3.561 MeV. On theoretical grounds, there are also self-consistent calcula¬ 
tions 1 ^ that suggest that Thomas-Ehrman shifts for nuclei near the drip line may be as 
large a few hundred keV. Because of the empirical nature of the INC interaction, however, it 
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is not clear how much of the Thomas-Ehrman effect has been absorbed into the interaction 
by £t. In addition, for nuclei near the drip-line, the estimate of the theoretical uncertainty 
is of the order of a 100-250 keV, and the effects of the Thomas-Ehrman shift are likely to 
he within the quoted uncertainties for the absolute binding energy. 

The parameters for the INC interaction used in this work are: €( 0 / 7 / 2 ) = 7.487 MeV, 
e(lp 3 / 2 ) = 7.312 MeV, €( 0 / 5 / 2 ) = 7.337 MeV, €( 1 ^ 1 / 2 ) = 7.240 MeV, Sc = 1.006, = 0.0, 

and = —4.2 x 10“^. Lastly, because of the difficulties associated with determining the 
0 / 5/2 and lp 3/2 single-particle energies, the uncertainties in the theoretical estimates of the 
^-coefficients for nuclei with A > 60 are increased from 33 keV to 45 keV. 


Shown in Table |ITl| are the results obtained for proton-rich nuclei whose binding energies 
are unknown in the mass region 46 < A < 69. The table lists the experimental binding 
energy of the neutron-rich analog, the predicted binding energy, one- and two-proton sepa¬ 
ration energies, as well as the Q-value for electron capture {Qec)- The ground-state spins 


were taken from Ref. and are also listed in the table. Wherever available, experimental 
binding energies tabulated in Ref. 0 were used in conjunction with theoretical values to 
compute Qec and the separation energies. 

In addition to the nuclei listed in Table B predictions for the “Coulomb” energy differ¬ 
ences for the the T = 1/2, A = 69 isodoublet and the T = 1, ^4 = 70 isotriplet, for which 
the binding energy of the neutron-rich member has not yet been measured experimentally. 


are given in Table |I^. The theoretical uncertainties for A = 70 include an uncertainty of 
20 keV in the c coefficient of Eq. (|l|). 

Shown in Fig. is a comparison of the binding energies reported here with those from 
three other theoretical studies. This comparison is illustrated by the difference ABE = 
i?£'(this work)—Ri7(other work), which is plotted in the hgure as a function of mass number, 
with the ordering the same as in Table B The error bars plotted at ABE = 0 represent 
the theoretical uncertainty of the binding energies listed in Table B The open circles show 
the comparison with the previous shell-model calculations of Ormand in Ref. 0 (A < 48), 
the open triangles the comparison with the binding energies of Cole in Ref. [§] {A < 52), 
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while the solid squares represent the comparison with the binding energies obtained from 
the unihed macroscopic-microscopic model of Moller and Nix |Q (46 < A < 70). While the 
results of Refs. are in agreement with those reported here, those of Moller and Nix are 
in severe disagreement for some nuclei. The origin of these differences is two fold. First, in 
the Moller-Nix study, the Coulomb-energy difference between analog nuclei is considerably 
smaller than in this work. This is illustrated in Fig. where the 6-coefficients for the nuclei 
listed in Table |T1| (open triangles) are plotted as a function of A in comparison with those 
derived from the Moller-Nix masses (solid squares). In the mass region 52 < ^4 < 64 the 
Moller-Nix 6-coefhcients are generally 100-200 keV smaller than those determined here. In 
addition, the Moller-Nix 6-coefficients are also in disagreement experimental trends, as is 
evidenced by the ®®Zn-^®Cu binding energy differences, where the Moller-Nix 6-coefhcient is 
9.683 MeV, which is 200 keV smaller than the experimental value of 9.881(40) MeV 0. The 
second reason for the large disagreement in Fig. can be attributed to poor reproduction 
of the mass of the neutron-rich analog nucleus. For example, the M511er-Nix ^®Cu mass 
excess is -54.8 MeV, which is in considerable disagreement with the experimental value of 
-56.3515(17) MeV tabulated in Ref. [0. 

Also shown in Fig. ^ is a comparison between the theoretical 6-coefficients and the 
systematic trends expected from the liquid-drop model of Eq. (P) (solid line) and the £t of 
Eq. (|I^) (dashed line). For the most part, the shell-model 6-coefhcients derived from the 
nuclei listed in Table H are in good agreement with the htted parameterization of Eq. (0. 
although they tend to be somewhat smaller than the systematic trend in the regions A < 50 
and A > 66. Note that for A < 50, this is a continuation of the trend for experimental data 
as is observed in Fig. |I] for 40 < A < 50. 

As a further illustration of the systematic trend for the Coulomb energy shifts, we exam¬ 
ine the halflives for the Fermi transition between analog J'^ = 0'^, T = 1 states in N = Z, 
odd-odd nuclei with A = 62, 66, and 70. The partial half-life for the (3 decay from the parent 
ground state to the state in daughter nucleus is given by 
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(13) 


where K = 2Tr^{\n2)h^/{mlc^) and K/Gy = 6170 ± 4 s |2^. The statistical rate function, 
fo^i, depends on the beta end-point energy and here it is evaluated using Eq. (10) of Ref. [0. 
For a pure Fermi transition between T = 1 analog states, the transition matrix element Mo^i 
is given by \J2{1 — 6 c), where 6 c is a small correction due to isospin-symmetry breaking. 
Recent calculations indicate that for these nuclei Sc is expected to be of the order 1-2%, 
and for the purpose of comparing with experimental data, will be taken to be equal to 1.5%. 
In general, Gamow-Teller transitions to excited states may also take place, and would tend 
to decrease the total beta-decay halflife. However, not only are the matrix elements for these 
transitions much smaller than for the Fermi transition, but since they occur to states in the 
daughter nucleus at a higher excitation energy, the statistical rate function is also much 
smaller. Hence, to a good approximation, the beta decay of these nuclei may be taken to be 
pure Fermi. Listed in Table 0 are the predicted Q-values for electron capture as well as a 
comparison between the experimental p4| , ^ , |30| and predicted beta-decay halflives. Given 
the fact that the statistical rate function strongly depends on the beta end-point energy (to 
the hfth power), the excellent agreement between the experimental and predicted beta-decay 
halflives for all three nuclei is a good indication that the overall systematic behavior of the 
Goulomb energy differences is well reproduced here. 

This section is concluded with a discussion on ®^As, which is important from an as- 
trophysical point of view. Because of the long beta-decay halflife for ®^Ge, if ®®As were 
signihcantly proton unbound, ®^Ge would then become a “waiting point” in the rp-process 
and would inhibit the production of heavier elements. If, however, the halflife of ®^As is 
dominated by beta decay, the rp-process will proceed through ®^As primarily by proton 
capture to ®®Se, although photodisintegration may begin to play an important role if ®®As 
is proton unbound . 


From Table |T|, ®®As is found in this work to be unbound to proton emission by 
0.428(254) MeV, with most of the uncertainty due to the uncertainty in the binding en- 
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ergy for (0.250 MeV). On the other hand, has been observed experimentally ||^ 
with a beta-decay halflife of lOOi'i^o® From the fact that no protons were observed in 
the stopping detector during this experiment, it may be inferred that the partial halflife 
for proton emission is signihcantly longer, and must be greater 1 s. The partial halflife for 
proton emission may be estimated using the WKB approximation, which is outlined in some 
detail in the next section [see Eq. and, in particular for proton emission, in Ref. 


A shell-model calculation within the /p-shell assuming a closed O/ 7/2 orbit and the FPD 6 * 
interaction yields 0.13 for the spectroscopic factor 6 *^. Using the potential parameters of 


Ref. 1 ^ 3 , a partial halflife for proton emission longer than 1 s requires the one-proton sep¬ 
aration energy to be greater than -0.23 MeV, which is in agreement with the value given 
in Table |T|. Because of the extreme sensitivity on the separation energy, however, it may 
never be possible to give a reasonable prediction for the partial halflife for proton emission 
without explicitly measuring the masses for both ®^Ge and ®^As. With the present uncer¬ 
tainty of 0.254 MeV, a range of 16 orders of magnitude is found for the halflife, i.e., between 
1.4 X 10“^^ to 1.6 X 10^ s. On the other hand, supposing that the binding energy of ®^Ge 
could be measured to within a few keV, a theoretical uncertainty of ~ 50 keV remains for 
®^As, which for a separation energy of -0.2 MeV, leads to a range of nearly four orders of 
magnitude in the proton partial halflife. 


IV. IDENTIFICATION OF CANDIDATES FOR DI-PROTON EMISSION 

In this section, the partial half-lives for di-proton emission are examined, with the inten¬ 
tion of identifying candidates amenable to experimental detection while taking into account 
theoretical uncertainties. As was mentioned in the introduction, the range of observable 
lifetimes for di-proton emission is limited by competing decay mechanisms and experimen¬ 
tal setups. In general, all candidates for di-proton emission have large /3-endpoints, and, as 
a consequence, the /3-decay halflives are expected to be of the order 1-100 ms []^. Also, in 
several experiments, such as in Ref. , the initial nucleus must live long enough to be iden- 
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tified. In this case, the limiting time is determined by the time-of-flight in the experimental 
apparatus. In general, these two conditions impose a practical limit on the observable halflife 


for di-proton emission to be in the range 10 ® — 10 


-3 


In Refs. p,|^, the di-proton decay halflives were estimated using r-matrix theory 
while taking the channel radius, Rq, to be 4 fm for all cases. In contrast, in Ref. the 
halflife for ^®Ni, was estimated using the Wentzel-Kramers-Brillouin (WKB) approximation. 
Because of uncertainties associated with the choice of the channel radius, the WKB approx¬ 


imation for the di-proton decay halflife will be used here. Following Ref. the WKB 
expression for the partial decay width is 


T 2 p = exp 

4p 


drk{r)) 


(14) 


where 6 “^ is the spectroscopic factor for hnding the di-proton in the correlated L = 0 state, /i 
is the reduced mass, rj„ and rout are the classical inner and outer turning points, respectively, 
the normalization factor Rf is determined by 




cos 


TT 


Uo 


dr'k{r') - 


= 1 , 


( 16 ) 


and k{r) is the wave number given by 




h m 


( 16 ) 


In Eq. ([T6|) , the asymptotic energy of the di-proton is Q 2 p = —*S' 2 p, V 2 p{r) is the average 
di-proton potential, and m*{r)/m is the proton effective mass. As in Ref. |^, V 2 p(r) is 
approximated by 2Vp(r), where Vp{r) is the self-consistent proton potential for the (Z — 
2, N) nucleus obtained from a Hartree-Fock or a Hartree-Fock-Bogoliubov calculation. Here, 
the halflives were computed using Hartree-Fock potentials using a Skyrme-type two-body 
interaction. It was found that the various Skyrme interactions give halflives that are in 
agreement to within a factor of two, and the results reported here were obtained using the 
Skyrme M* interaction. In addition, the halflives computed using Eq. (0 were found to 
be approximately an order of magnitude shorter than those obtained using the r-matrix 
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representation with i?o = 4 fm (as was used in Refs. ^,0). On the other hand, if the 
channel radius is chosen to be equal to the classical inner turning point, Tin, the r-matrix 
approach yields halflives that are within a factor of two of the WKB method. 

The spectroscopic factor, 6 , can be evaluated within the framework of the shell model. 
For di-proton emission the spectroscopic factor can be estimated using the cluster-overlap 
approximation |^, namely 




(17) 


where fc. A, and are parameters dependent on the model space and the emitted cluster, 
and 'ipc is a two-proton cluster wave function in which the relative motion of the particles is 
governed by the OS' state, and is obtained by diagonalizing an SU3 conserving interaction 
within the shell-model conhguration space . 


Of all the quantities in Eq. (|^), the di-proton decay rate is most sensitive to the two- 
proton separation energy S 2 p. Indeed, it was illustrated in Ref. 0] that an uncertainty of 
±100 keV in a separation energy of the order 500 keV can lead to a range of nearly six 
orders of magnitude in the di-proton decay halflife. In contrast, the spectroscopic factors 
are expected to be of the order 0.5-0.75 |^, and shouldn’t lead to any more than an order 
of magnitude decrease in the decay rate (increase in the halflife). Given that the theoretical 
uncertainties in the separation energy for each of the di-proton emitters considered in this 
work are all greater than 175 keV, an accurate estimate of the spectroscopic factor is not 
needed in order to obtain an order-of-magnitude estimate of the di-proton halflife for the 
purpose of identifying the best candidates for experimental observation. Hence, the lifetimes 
reported here are evaluated assuming 9“^ = 1 with the understanding that they are probably 
too short by a factor of two to four. 

Listed in Table Q are the halflives (ti /2 = hln 2 /r 2 p) associated with di-proton emission 
for all nuclei in Table H that are predicted to be unstable to two-proton emission while 
being bound to proton emission. Also, for the purpose of comparison, the halflives for ^®Ti, 
^^Fe, and ^®Ni given in Ref. 0 are also listed in the table. Given the practical limitations on 
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the halflife for the experimental observation of this decay mode, the best candidates are ‘^^Fe, 
^®Ni, and ®^Se. Of these three, perhaps the best is ^^Fe since it is likely that it has already 
been identihed experimentally [Q. On the other hand, both ®®Ge and ®^Kr have halflives 
that are long enough to make them marginal candidates for experimental observation. 

V. CONCLUSIONS 

In this work. Coulomb energy differences between mirror nuclei with 46 < A < 70 
were computed within the framework of the nuclear shell model using an effective Coulomb 
plus isotensor interaction. Absolute binding energies for proton-rich nuclei are predicted by 
adding the Coulomb energy differences to the experimental binding energy of the neutron- 
rich analog. With these binding energies, proton separation energies are computed, and the 
location of the proton drip-line is delineated. 

The computed Coulomb energy differences were also compared with systematic trends 
predicted by the liquid-drop model and a £t to experimental 6-coefficients assuming a 
dependence. It was found that the shell-model calculations were in good agreement with the 
systematic trends, except for A < 50 and A > 66. As a further test on the systematic trend 
of the shell-model Coulomb shifts, halflives for the Fermi transitions in odd-odd, N = Z 
nuclei with A = 62, 66, and 70 were computed and found to be in excellent agreement with 
experimental data. The shell-model binding energies predicted here were also compared with 
three previous works. While the results of Ormand (only for A < 48) and Cole [§| (only 
for A < 52) are in good agreement with those reported here, those of Moller and Nix are 
not. It was found that the disagreement with the Moller-Nix masses is due to differences 
in both the Coulomb energy shifts and the binding energy of the neutron-rich analog. For 
the most part, the data presented in Fig. is the only data that is explicitly sensitive to 
a parameterization of the Coulomb energy. Civen the importance of analog symmetry and 
the overall success of the IMME, any global parameterization of binding energies should 
include a proper description of the Coulomb energy differences. Towards this end, perhaps 
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the best approach is to determine the parameters of a microscopic-macroscopic model using 
the neutron-rich binding energies, while hxing the parameters of the Coulomb plus isovector 
part so as to reproduce the Coulomb energy shifts between mirror nuclei. Even in this limit, 
however, it has to be noted that the systematic parameterization is capable of reproducing 
the experimental b coefficients of the IMME only at the level of approximately 100 keV. 

Finally, two-proton separation energies were also computed, and halflives associated with 
correlated di-proton emission were computed using the WKB approximation. Given practi¬ 
cal constraints on the halflife for the observation of this decay mode imposed by competition 
with beta decay and experimental setups, the best candidates for experimental observation 
are predicted to be ^^Fe, "^®Ni, and ®^Se. 
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TABLES 



TABLE 1. 

Values of hto used for /p-shell nuclei 


A 

huj (MeV) 

A 

huj (MeV) 

40 

10.603 

60 

10.156 

41 

10.603 

61 

10.087 

42 

10.603 

62 

10.017 

43 

10.608 

63 

9.954 

44 

10.614 

64 

9.890 

45 

10.603 

65 

9.786 

46 

10.592 

66 

9.681 

47 

10.581 

67 

9.589 

48 

10.570 

68 

9.496 

49 

10.560 

69 

9.460 

50 

10.550 

70 

9.424 

51 

10.539 

72 

9.331 

52 

10.528 

73 

9.168 

53 

10.507 

74 

9.203 

54 

10.486 

76 

9.032 

55 

10.470 

77 

9.100 

56 

10.454 

78 

8.923 

57 

10.376 

79 

9.869 

58 

10.298 

80 

8.816 

59 

10.227 
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TABLE II. Comparison between theoretical (with FPD6* and experimental excitation energies 


(in MeV) of the first = 2~^ state in even-even N = Z fp-shell nuclei 



Expt. 

EPD6* 

60 Zn 

1.004“ 

0.825 

o^Ge 

0.902^ 

0.700 

6®Se 

0.854^ 

0.600 

72Kr 

0.709^ 

0.707 

76sr 

0.26R 

0.752 

80 zr 

0.289^ 

- 


“from Ref. 
^from ref. 


2 § 

|26|1 
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TABLE III. Predicted binding energies, one- and two-proton separation energies {Sp and S 2 p, 
respectively), /3-decay end-point energies for proton rich nuclei with 46 < A < 70. The absolute 
binding energies were computed with theoretical Coulomb energy shifts added onto the experimen- 


tal binding energy for the neutron-rich analog, also listed 

in the table. 






BEthy 


TD jpanaloq 
^^exp 

Sp 

S2p 

Qec 


T, 


(MeV) 

"^Z-analog 

(MeV) 

(MeV) 

(MeV) 

(MeV) 

46 Mn 

2 

4+ 

364.186(132) 

Sc 

396.610(1) 

0.156(146) 

3.234(139) 

17.007(134) 

46Fe 

3 

0+ 

350.144(198) 

46Ca 

398.769(2) 

1.408(224) 

0.328(215) 

13.260(238) 

^^Mn 

3/2 

5/2- 

382.326(99) 

47Ti 

407.072(1) 

0.351(101) 

5.237(100) 

12.020(100) 

^■^Fe 

5/2 

7/2- 

365.973(165) 

^"Sc 

407.254(2) 

1.787(211) 

1.943(177) 

15.571(192) 

47Co 

7/2 

7/2- 

348.349(231) 

47Ca 

406.045(2) 

-1.795(304) 

-0.387(254) 

- 

^^Mn 

1 

4+ 

397.101(66) 

48 V 

413.904(3) 

1.973(68) 

6.740(66) 

13.579(66) 

00 

2 

0+ 

385.106(132) 

48Ti 

418.698(1) 

2.780(165) 

3.131(134) 

11.213(148) 

00 

o 

O 

3 

6+ 

365.153(198) 

Sc 

415.487(5) 

-0.820(258) 

0.967(238) 

- 

48Ni 

4 

0+ 

348.854(264) 

48Ca 

415.991(4) 

0.505(351) 

-1.290(330) 

15.517(330) 

49Fe 

3/2 

7/2- 

399.802(99) 

49v 

425.457(1) 

2.701(119) 

4.674(100) 

12.963(102) 

^^Co 

5/2 

7/2- 

384.184(165) 

49Ti 

426.841(1) 

-0.922(211) 

1.858(192) 

- 

49 Ni 

7/2 

7/2- 

365.830(231) 

49 Sc 

425.618(4) 

0.677(304) 

-0.143(284) 

17.572(284) 

49Cu 

9/2 

3/2- 

344.413(297) 

49Ca 

421.138(4) 

-4.441(397) 

-3.936(376) 

- 

5°Co 

2 

6+ 

400.060(132) 

50v 

434.790(1) 

0.258(165) 

2.959(148) 

16.585(145) 

50Ni 

3 

0+ 

385.693(198) 

50Ti 

437.780(1) 

1.509(258) 

0.587(238) 

13.585(238) 

50Cu 

4 

5+ 

362.299(264) 

“Sc 

431.674(16) 

-3.531(351) 

-2.854(330) 

- 

50Zn 

5 

0+ 

340.823(330) 

“Ca 

427.491(9) 

-3.590(444) 

-8.031(423) 

- 

54Co 

3/2 

7/2- 

417.864(99) 

5icr 

444.306(1) 

0.164(116) 

4.317(102) 

12.868(100) 

51 Ni 

5/2 

7/2- 

401.684(165) 

Sly 

445.841(1) 

1.624(211) 

1.882(192) 

15.398(192) 

s^Cu 

7/2 

3/2- 

382.472(231) 

“Ti 

444.153(1) 

-3.221(304) 

-1.712(284) 

- 

“Co 

1 

6+ 

432.912(66) 

52 Mn 

450.851(2) 

1.398(68) 

6.283(66) 

14.003(67) 


25 






52 Ni 

2 

0+ 

420.478(132) 

52Gr 

456.345(1) 

2.614(165) 2.778(145) 

11.652(148) 

52CU 

3 

3+ 

399.399(198) 

52y 

453.152(1) 

-2.285(258) -0.661(238) 

- 

52Zn 

4 

0+ 

380.321(264) 

52Ti 

451.961(7) 

-2.151(351) -5.372(330) 

- 

53 Ni 

3/2 

7/2- 

435.558(99) 

53 Mn 

462.905(2) 

2.646(119) 4.044(100) 

12.956(101) 

53CU 

5/2 

3/2- 

418.835(165) 

53Gr 

464.285(2) 

-1.643(211) 0.971(192) 

- 

53 Zn 

7/2 

7/2- 

397.948(231) 

SSy 

461.631(3) 

-1.451(304) -3.736(284) 

- 

54Cu 

2 

3+ 

434.906(132) 

54Mn 

471.844(2) 

-0.652(165) 1.994(148) 

- 

54 Zn 

3 

0+ 

418.605(198) 

54 Gr 

474.004(1) 

-0.230(258) -1.873(238) 

15.519(238) 

54Ga 

4 

3+ 

393.891(264) 

54v 

467.744(15) 

-4.057(351) -5.508(330) 

- 

55Cu 

3/2 

3/2- 

452.997(99) 

55Fe 

481.057(1) 

-0.153(111) 3.701(101) 

13.568(100) 

55 Zn 

5/2 

5/2- 

435.071(165) 

55 Mn 

482.071(1) 

0.165(211) -0.487(192) 

17.144(192) 

55Ga 

7/2 

3/2- 

414.644(231) 

55 Gr 

480.250(1) 

-3.961(304) -4.191(284) 

- 

56Cu 

1 

4+ 

467.899(66) 

56Co 

486.906(1) 

0.552(67) 

5.166(66) 

15.307(67) 

56 Zn 

2 

0+ 

454.214(132) 

sepe 

492.254(1) 

1.217(165) 1.064(141) 

12.903(148) 

56Ga 

3 

3+ 

432.226(198) 

56 Mn 

489.315(1) 

-2.845(258) -2.680(238) 

- 

56Ge 

4 

0+ 

412.381(264) 

56Gr 

488.507(10) 

-2.263(351) -6.224(330) 

- 

57Zn 

3/2 

7/2- 

469.440(99) 

57Co 

498.282(1) 

1.541(119) 2.093(100) 

14.461(100) 

57Ga 

5/2 

3/2- 

451.874(165) 

57Fe 

499.885(1) 

-2.340(211) -1.123(192) 

- 

57Ge 

7/2 

5/2- 

430.634(231) 

57Mn 

497.992(3) 

-1.592(304) -4.437(284) 

- 

O 

00 

2 

2+ 

468.039(132) 

58Co 

506.855(2) 

-1.401(165) 0.140(148) 

- 

58Ge 

3 

0+ 

451.578(198) 

58Fe 

509.945(1) 

-0.296(258) -2.636(238) 

15.679(238) 

58 As 

4 

3+ 

426.697(266) 

58 Mn 

504.480(30) 

-3.937(352) -5.529(331) 

- 

59Ga 

3/2 

3/2- 

486.040(99) 

59Ni 

515.453(1) 

-0.920(111) 1.357(100) 

- 

59Ge 

5/2 

7/2- 

468.097(165) 

59Co 

517.308(1) 

0.058(211) -1.343(192) 

17.161(192) 

59 As 

7/2 

3/2- 

447.648(231) 

59Fe 

516.526(1) 

-3.930(304) -4.226(284) 

- 

69Ga 

1 

2+ 

500.080(66) 

59Cu 

519.933(3) 

0.080(77) 

2.971(66) 

14.130(67) 

69Ge 

2 

0+ 

487.127(132) 

60Ni 

526.842(1) 

1.087(165) 0.167(141) 

12.171(148) 
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6° As 

3 

5+ 465.094(198) 

60Co 

524.800(1) -3.003(258) -2.945(238) 

- 

6iGa 

1/2 

3/2- 515.179(48) 

64 Zn 

525.223(16) 0.187(49) 5.307(48) 

9.262(50) 

Ge 

3/2 

3/2- 501.415(135) 

64Cu 

531.642(2) 1.335(150) 1.415(141) 

12.982(143) 

61 As 

5/2 

3/2- 484.381(225) 

6lNi 

534.595(1) -2.746(261) -1.659(246) 

- 

62Ge 

1 

0+ 517.720(91) 

62 Zn 

538.119(10) 2.541(102) 2.728(91) 

9.664(95) 

62 As 

2 

1+ 499.816(180) 

62CU 

540.529(4) -1.599(225) -0.264(192) 

- 

62Se 

3 

0+ 484.239(270) 

62Ni 

545.259(1) -0.142(351) -2.888(301) 

14.795(325) 

63Ge 

1/2 

3/2- 530.597(110) 

63Ga 

540.930(100) 2.431(113) 5.374(111) 

9.551(148) 

63 As 

3/2 

3/2- 516.321(135) 

63 Zn 

547.232(2) -1.399(163) 1.142(143) 

- 

63Se 

5/2 

3/2- 499.885(225) 

63Cu 

551.382(1) 0.069(288) -1.530(262) 

15.654(262) 

64 As 

1 

0+ 530.315(90) 

64 Ga 

551.147(4) -0.282(142) 2.149(94) 

14.853(266) 

64Se 

2 

0+ 517.411(180) 

64 Zn 

559.094(2) 1.090(225) -0.309(202) 

12.122(201) 

65 As 

1/2 

3/2- 545.522(46) 

65Ge 

556.010(10) -0.428(254) 4.592(110) 

- 

65Se 

3/2 

3/2- 531.473(135) 

66Ga 

563.036(2) 1.158(162) 0.876(174) 

13.267(143) 

65 Br 

5/2 

1/2- 514.580(225) 

66 Zn 

567.020(2) -2.831(288) -1.741(262) 

- 

66Se 

1 

0+ 548.091(95) 

66Ge 

569.290(30) 2.569(105) 2.141(267) 

10.087(112) 

66 Br 

2 

0+ 529.780(180) 

66Ga 

572.176(3) -1.693(225) -0.535(201) 

- 

66Kr 

3 

0+ 514.579(270) 

66Zn 

578.133(2) -0.001(351) -2.832(325) 

14.419(325) 

67Se 

1/2 

5/2- 560.882(110) 

67 As 

571.610(100) 1.922(125) 4.872(110) 

- 

67Br 

3/2 

1/2- 546.355(135) 

67Ge 

578.398(5) -1.736(165) 0.833(143) 

- 

67Kr 

5/2 

1/2- 529.935(225) 

67Ga 

583.403(2) 0.155(288) -1.538(262) 

15.638(262) 

68 Br 

1 

3+ 560.365(135) 

68 As 

581.910(100) -0.517(174) 1.405(147) 

- 

68Kr 

2 

0+ 547.668(180) 

68Ge 

590.792(6) 1.313(225) -0.423(204) 

11.915(225) 

68Rb 

3 

1+ 526.980(270) 

68Ga 

591.680(2) -2.955(351) -2.800(325) 

- 

69 Br 

1/2 

3/2- 575.737(54) 

69Se 

586.620(30) -0.663(305) 4.127(114) 

- 

69Kr 

3/2 

5/2- 561.477(138) 

6® As 

594.180(30) 0.075(193) -0.442(176) 

14.515(148) 

™Rb 

2 

4+ 559.398(187) 

76 As 

603.520(50) -1.042(232) -0.967(230) 

- 
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TABLE IV. Predictions for the “Coulomb” energy difference for the T = l/2, A = 69 and 


T = 1, A = 70 nuclei. 




ABE (MeV) 

®9Br-69Se 

3/2- 

-10.883(45) 

70Kr-70Br 

0+ 

-11.241(50) 

70Br-70Se 

0+ 

-10.801(50) 


TABLE V. Comparison between experimental and predicted beta-decay halflives for odd-odd, 
N = Z Eermi transitions. The predicted (5 EC-value is also given. 



Qec (MeV) 

ti/2 (ms) 

tT/l (ms) 

62Ga 

9.191(50) 

115(2) 

116.1(2)“ 

66 As 

9.592(50) 

94(2) 

95.8(2)^ 

70Br 

10.019(50) 

76(3) 

79.1(8)“ 


From Ref. 

29 

From Ref. 


From Ref. 
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TABLE VI. Half-lives for di-proton emitter candidates. Also listed are the theoretical predic¬ 


tions for the one- and two-proton separation energies 



Sp (MeV) 

S 2 p (MeV) 

H/2 (s) 

tr/T (s) 

+max 

^1/2 

Ref. [7] 

38 Ti 

0.438(164) 

-2.432(132) 

9 X 10"36 

4 X 10-36 

2 X 10-36 

45 Fe 

-0.010(198) 

-1.279(181) 

10-6 

10-3 

10-4 

48Ni 

0.502(164) 

-1.137(210) 

3 X 10-3 

10-6 

4 

This work 

48Ni 

0.505(351) 

-1.290(330) 

4 X 10-6 

5 X 10-9 

0.09 

49Ni 

0.677(304) 

-0.143(284) 

3 X 10^® 

6 X 1034 

oo 

55Zn 

0.165(211) 

-0.487(192) 

5 X 1034 

2 X 109 

3 X 1030 

59Ge 

0.058(211) 

-1.343(192) 

10-3 

10-6 

0.3 

63 Se 

0.069(288) 

-1.530(262) 

6 X 10-6 

3 X lO-"^ 

5 X 10-2 

64Se 

1.090(225) 

-0.309(202) 

5 X 1032 

6 X 1033^ 

4 X 103’9 

66Kr 

-0.001(351) 

-2.832(325) 

3 X 10-32 

2 X 10-33 

6 X 10-33 

67Kr 

0.155(288) 

-1.538(262) 

2 X 10-3 

10-6 

0.2 

68Kr 

1.313(225) 

-0.423(204) 

3 X 1024 

8 X 1033 

5 X 1049 

69Kr 

0.075(193) 

-0.442(176) 

2 X 1023 

2 X 1034 

1040 
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FIGURES 

FIG. 1. Dependence of 6-coefficients as a function of mass number A. Experimental data are 
represented by the solid squares, while the values from the liquid-drop formula and the fit [Eqs. ® 
and (0)] are represented by the solid and dashed lines, respectively. 


FIG. 2. Difference between absolute binding energies listed in Table 111 with those of Moller 


and Nix [5] (solid squares), Ormand [7], and Cole [8] (open triangles). The error bars at ABE = 0 


denote the theoretical uncertainty of the binding energies listed in Table 111. 


FIG. 3. Dependence of theoretical 6-coefficients for nuclei listed in Table (open triangles) 
as a function of mass number A. For comparison, the 6-coefficients derived from the unified mi¬ 
croscopic-macroscopic model of Moller and Nix [5] are also shown (solid squares). The systematic 


behavior as expected from Eqs. (|^ and (|lb[) are represented by the solid and dashed lines, respec¬ 
tively. 
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T 



A 


—■— Moller-Nix [5] 
o Ormand [7] 

^ Cole [8] 
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